Turbulence and coarsening in active and passive binary mixtures 
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Phase separation between two fluids in two- dimensions is investigated by means of Direct Nu- 
merical Simulations (DNS) of coupled Navier-Stokes and Cahn-Hilliard equations. We study the 
phase ordering process in the presence of an external stirring acting on the velocity field. For both 
active and passive mixtures we find that, for a sufficiently strong stirring, coarsening is arrested in a 
stationary dynamical state characterized by a continuous rupture and formation of finite domains. 
Coarsening arrest is shown to be independent of the chaotic or regular nature of the flow. 



PACS numbers: 47.27.-i,05.70.Ln 

When a binary fluid mixture at the critical concentra- 
tion is cooled from a high temperature to a sufficiently 
low temperature (below a critical one), the original ho- 
mogeneous phase becomes unstable and spontaneously 
evolves into two phases separated by an interface. As 
time advances an out-of-equilibrium process of phase or- 
dering takes place through the formation of domains of a 
single phase that grow algebraically in time [1]. In fluids, 
the presence of a hydrodynamic velocity field makes this 
process more complicated than the corresponding one in 
solid alloys. For instance, since Siggia's seminal work [2], 
it is well known that hydrodynamics may accelerate the 
domain growth (see Refs. [3] and [4] for recent develop- 
ment in three and two dimensional fluids, respectively). 
Phase ordering dynamics becomes even more complex 
and less understood when the fluid mixture is externally 
driven [5, 6]; beyond their theoretical interest, phase sep- 
arating binary fluids under flow embody a great techno- 
logical interest [7] for their distinctive rheological prop- 
erties. This problem has been extensively investigated 
in shear flows [8, 9, 10, 11] where coarsening becomes 
highly anisotropic: the single-phase domain growth ac- 
celerates in the shear direction while in the transversal 
one the growth is arrested [10, 11] or strongly slowed 
down [9]. Less clear is the case in which the mixture is 
stirred by a turbulent flow [7, 12, 13, 14]. Here, phase 
separation may be completely suppressed [13], or a dy- 
namical steady state with domains of finite length and 
well defined phases may develop [6, 14]. A similar phe- 
nomenology has been experimentally observed in stirred 
immiscible fluids [15]. This Letter aims to clarify the na- 
ture of the non-equilibrium steady state, characterized 
by the continuous rupture and formation of domains. 

Previous investigations focused on passive binary mix- 
tures (when the feedback of the phase ordering on fluid 
velocity is neglected) in random flows [16] and in chaotic 
flows [17] (in a Lagrangian sense, i.e. two initially very 
close particles separate exponentially in time [18]). Here 
we focus on the phase ordering dynamics of active two- 



dimensional binary mixtures in which the fluid is driven 
by chemical potential inhomogeneities [6]. By means of 
numerical experiments we show that coarsening arrest 
is a generic and robust phenomenon, whose existence 
can be understood by an energy conservation argument. 
Moreover, we show that in the passive limit Lagrangian 
Chaos is not necessary for coarsening arrest. 

In the presence of stirring, the main question concerns 
the competition between thermodynamic forces, driving 
the phase segregation, and fluid motion, leading to mix- 
ing and domains break-up. For very high flow intensities 
phase separation can be completely suppressed [12, 13] 
due to mixing of the components and inhibition of in- 
terface formation. In active mixtures with very low vis- 
cosities such a phenomenon may be self-induced by the 
feedback [4, 19]: the fluid responds vigorously to local 
chemical potential variations and re-mixes the compo- 
nents. On the other hand, stirring may lower the critical 
temperature [5, 14, 20]. However, by performing a deeper 
quench, phase separation in a nontrivial statistically sta- 
tionary state may still develop [6]. 

Being interested in deep quenching, here we work at 
zero temperature, as in Refs. [16, 17]. We consider a 
symmetric (50% — 50%) mixture of two incompressible 
fluids of equal density p = 1 and viscosity u. Such a 
bi-component system is described by a scalar order pa- 
rameter 0{r^t)^ the local fraction of the two fluids. The 
associated Landau-Ginzburg free energy reads [1]: 

^0]=Jdr(^-^e' + \e' + ^-^\v0\''^ , (1) 

being ^ the equilibrium correlation length, which pro- 
vides a measure of the interface width. The dynamics is 
then governed by the Cahn-Hilliard (CH) equation 

S,^ + i; • = FV" ^ = F V V (2) 

where /i = —0-\-0^ — ^^V^^ is the chemical potential, and 
F is a mobility coefficient that we assume constant and in- 
dependent of 0. Hydrodynamics enters in Eq. (2) through 
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the convective term. The order parameter is transported 
by the two-dimensional velocity field which evolves 
according to the Navier-Stokes (NS) equation 



dtv -Vv = vV'^v -Vp- 6Vfi + / . 



(3) 



where p is the pressure. The fiuid is forced by the ex- 
ternal mechanical force / and by local chemical poten- 
tial variations — ^V/i (the two phases want to demix and 
thus force the fiuid, see also Ref. [3] for a detailed deriva- 
tion). This latter term can be rewritten as — ^^V^^V^ 
plus a gradient term which can be absorbed into the pres- 
sure [6]. Therefore Eq.(3) formally reduces to the 2D 
magnetohydrodynamics (MHD) equation for the veloc- 
ity field. Actually phase ordering and MHD share many 
phenomenological properties [6] . 

We numerically integrate the coupled equations (2)- 
(3) by means of a standard pseudo-spectral code imple- 
mented on a two-dimensional periodic box of size 27r x 27r 
with 512^ collocation points. Statistical analysis of do- 
main sizes is obtained by considering the characteristic 
length, defined as = ((1 — 6>^))~^ where (...) denotes 
spatial average [21]. The initial condition for the order 
parameter is a high temperature configuration with set 
as white noise in space. In the presented results time is 
rescaled with the diffusive time tm = ^^/r. 

Unstirred case. Starting from the initial configuration 
with the fiuid at rest (v = 0), after a few diffusive time 
scales tm, sharp interfaces appear and phase separation 
proceeds through domain coarsening. At long times, the 
domains length L - the only characteristic scale of the 
system (provided L - grows in time as a power 

law [1] . In 2D different regimes have been predicted and 
observed [4]: L{t) ~ t^/^, as in fiuids at rest, for high 
viscosity; L{t) ~ t^/^ for lower viscosities. At intermedi- 
ate values of viscosity it is still unclear whether there is 
only one characteristic scale [4]; for <C 1 and low mo- 
bility r <C 1 mixing may overwhelm phase demixing [19]. 
In the following we will limit our analysis to the turbu- 
lent, low viscous, regime where scaling exponent 2/3 is 
expected. This exponent can be dimensionally derived 
by balancing the inertial term v ■ with ^V/i in (3), 
and assuming that L{t) is the only length scale of the 
system. The scaling behavior of L{t) implies the follow- 
ing ones for the kinetic energy and the enstrophy [3]: 
/C = {v^)/2 - and Q = {u;^)/2 - t'^/^ {uj = V x v 

is the vorticity). Fig. 1 shows that the scaling predictions 
are well reproduced by our DNS. We remark that in the 
absence of stirring, phase separation is accelerated by the 
presence of hydrodynamics. 

Stirred case. We now consider the presence of an 
external mechanical forcing acting on the velocity field. 
As customary in turbulent simulations, energy is injected 
by means of a random, time uncorrelated, homogeneous, 
and isotropic Gaussian process with amplitude F which 
is restricted to a few Fourier modes around kf (this iden- 
tifies the injection scale if ~ 27r//c/). The ^-correlation 
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FIG. 1: (color online) (a) L{t) vs t obtained by DNS of (2) 
and (3) with ^ = 0.015, u — 10~^, without external forcing, 
(b) Kinetic energy AC = ('i^^)/2 (bottom) and enstrophy = 
(a;^)/2 (top) vs t in the same run. 



in time allows for controlling the kinetic energy input 
Cin = F'^Uf (being n/ the number of excited Fourier 
modes). Equations (2)- (3) are integrated starting with 
i; = 0. In Fig. 2 we show typical snapshots of the order 
parameter at different fiow intensities, and for two ran- 
dom forcing with ^/ ~ 26^ and if ^ 84^. As it is evident, 
the stronger the stirring intensity the smaller the typical 
domain length. This is confirmed by the temporal evolu- 
tion of L(t) at varying the external forcing (Fig. 3). After 
an initial growth characterized by the 2/3 scaling expo- 
nent, L(t) stabilizes at a value L* that decreases with the 
stirring intensity. Both the kinetic energy ]C(t) and the 
enstrophy Q.(t) (not shown here) stabilize at correspond- 
ing values. Therefore a well defined statistically steady 
state is reached. 

A closer inspection of Fig. 2 reveals qualitative dif- 
ferences in domain shapes. When L* is larger than the 
forcing scale if (left) the domains are almost isotropic, 
while in the case L"" < if (right) the underlying veloc- 
ity field reveals itself through the filamental structure of 
the domains. Nevertheless, coarsening process is always 
arrested confirming the robustness of the phenomenon. 
Stirring always selects a scale through the competition 
between the thermodynamic forces and the stretching in- 
duced by local shears that deform and break the domains. 
Estimating the shear rate as 7 = Urms/L'' we find that. 





FIG. 2: Snapshots of at time t — 4000 at varying the forc- 
ing intensity F with = 26^ (left) and ^/ = 84^ (right). 
Black/white codes — ^1. Other parameters are set as in 
Fig. 1 
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FIG. 3: (color online) L vs t at varying F, from top F — 
(thick curve) and F=0.05, 0.10, 0.15, . . . , 0.30. Data refer to 
DNS with if = 84:^ (the case with £f ^ 26^ is qualitatively 
similar). The straight line displays the scaling t^^^ . Inset: L* 
vs Urms/L*'', the straight line has slope —0.29, point size is of 
the order of the statistical error. 
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FIG. 4: (color online) Results of DNS in the passive case, (a) 
L vs t with (from top) f5 — 0.25,0.5, 1,2,4. Forcing param- 
eters are kf — 70 and F = 3 x 10~^. The straight line has 
slope 1/3. (b) Same as (a), in linear scale, for a frozen velocity 
field (see text), from top /3 = 1,2,4, 10. Inset: L* vs the 
straight line corresponds to L* ~ ^ point size is of the 

order of the statistical error. DNS were performed with hyper 
dissipation —V2^v with V2 — 10~^, and friction coefficient 
a = 0.1. The parameters of Eq. (2) are ^ = 0.015, F ^ 0.02. 



for the case with L* < L* ~ ^-0.29 (^^^^ m_^Q\, of 
Fig. 3), in fairly good agreement with experiments and 
simulations in pure shear flows [11]. However we should 
mention that in our settings, homogeneous and isotropic 
flows, there is not a well defined rate 7 as in genuine 
shear flows. The definition adopted here is a dimensional 
estimation of the shear rate at the arrest scale. In the 
case L* > no clear scaling behavior is observed. 

The existence of a stationary state can be understood 
in terms of conservation laws. Due to the presence of two 
inviscid quadratic invariants, /C and 1], the single fluid 2D 
NS equation (i.e. (3) without the feedback term) is char- 
acterized by a double cascade [22]: /C flows toward the 
large scales {r>if) and Q toward the small ones {r<if). 
By switching on the coupling term, the following balance 
equation for the total energy E = JC -\- ^ holds [23] 

^=-^(|Vt;|')-r(|VA.p)+e,„. (4) 

Now it is clear that in the unstirred case (e^n =0) asymp- 
totically an equilibrium state, corresponding to fluid 
at rest v = and complete phase separation (minimum 
of free energy), is reached. In this case the velocity 
has only a transient role, determining the scaling of 
the coarsening process. On the contrary if e^^ 7^ a 
nontrivial stationary state stems from the balance of 
dissipative and input terms in the r.h.s. of (4). 

It is worth mentioning that when the stirring intensity 
becomes high enough to overcome the feedback term, the 
kinetic energy dissipation induced by the |V/ip term is 
no more effective. Indeed, when ^/ is much larger than ^ 
and F is very high, the coupling term becomes negligible 
and we observe the single-fluid phenomenology with an 
inverse energy cascade. 

Passive binary mixtures We now consider the case in 
which the coupling term in (3) is switched off and conse- 
quently the order parameter is passively transported by 
the velocity field. This case has been already considered 



in [16, 17]. In order to obtain a statistically stationary 
state, as customary we added a large scale friction term 
—av to the Navier-Stokes equation [24]. The velocity 
field in (2) is rescaled by a factor /3; this is a numerically 
convenient way to change the velocity intensity and to 
study the effect of stirring on coarsening. For /3 = 
Eq. (2) recovers the Cahn-Hilliard equation in a fluid at 
rest for which L{t) ~ t-^^^. For /3 > we observe the 
following phenomenology (Fig. 4a). For small values of 
/3 (weak stirring) we did not find a clear evidence of 
coarsening arrest. This is likely due to finite size effects 
hiding the phenomenon, i.e. L* becomes comparable 
with, or even larger than, the box size. For /3 large 
enough (strong stirring), the existence of an arrest scale 
L* (that decreases with /3) is well evident. 
Previous studies stressed the importance of Lagrangian 
Chaos in the coarsening arrest phenomenon [17]. 
Now, in order to elucidate this point, we discuss a 
non-chaotic example. It is well known (see, e.g., [18]) 
that two-dimensional stationary flows do not generate 
chaotic trajectories. We have thus integrated (2) in 
a frozen configuration of the turbulent velocity field: 
vix^t) = v{x). As shown in Fig. 4b, domain growth 
is strongly weakened and finally arrested, even in this 
non-chaotic flow. For moderate velocity intensities, 
L{t) still grows in time, but with a much slower scaling 
law than the dimensional prediction for fluids at rest, 
t^/^. This slowing down is probably due to a different 
growth mechanism: after an initial transient, a slow 
process of droplet passage among close domains is 
indeed observed. However, for high enough intensities a 
complete stabilization of the domain length is realized. 
This suggests that the main ingredient for coarsening 
arrest is the presence of local shears that overwhelm the 
surface tension driving force. The dependence of L* on 
the shear rate (here naturally defined as P) is shown in 
the inset of Fig. 4b. We find a power law behavior with 
exponent —0.28 very close to the one observed in the 
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FIG. 5: (left) Snapshots of at t=6xl0^ for the cellular flow 
with different C/ and K = 8, here ^ = 0.018, r = 0.1. (right) 
L vs t, from top to bottom C/ = 0, 0.125, 0.25, 0.5, 1.0, 2.0, 4.0. 
Inset: the same in linear scale. 

active case and in shear flows [11], while in the chaotic 
case no clear scaling is observed. 

To further support the marginal role of La- 
grangian Chaos in coarsening arrest, we report the 
results obtained in a stationary regular cellular flow 
{vx = sin(Kx) cos(K7/) , Vy = —U cos{Kx) sm{Ky)) 
where U fixes the velocity amplitude and K the char- 
acteristic scale. As shown in Fig. 5 (left), for large 
intensities the order parameter is frozen into a random 
chessboard pattern with a finite length. At lower intensi- 
ties a growth much slower than in the absence of the flow 
is still visible (Fig. 5 (right)) coming from a slow droplet 
migration from one cell to another. At large [/'s the 
shear between the counter-rotating vortexes overwhelms 
the demixing induced by the thermodynamic forces, 
breaking the domains which freeze into the cells. 

In conclusion, we have shown that in the presence 
of an external stirring coarsening process is slowed 
down both for active and passive mixtures. We have 
also demonstrated that the phenomenon of coarsening 
arrest, first predicted in [6], does not necessarily require 
a chaotic fiow, as suggested in [17], but is a consequence 
of the competition between thermodynamic forces and 
stretching induced by local shears. Our investigation in 
both active and passive mixtures shows that this be- 
havior is robust. Moreover we found numerical evidence 
that the dependence of the arrest scale on the shear 
rate follows a power law behavior with an exponent 
close to the one measured in experiments and numerical 
simulations in pure shear fiows [11]. Our results might 
suggest the existence of a mechanism independent of 
the nature of the flow in the coarsening arrest. Further 
numerical and experimental investigations, with the aim 
of clarifying the dependence of the arrest scale on the 
flow properties, would be extremely interesting. 
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